Numerical simulation of the twist-grain-boundary phase of chiral liquid crystals 
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We study the structure of the twist-grain-boundary phase of chiral liquid crystals by numerically 
minimizing the Landau-de Gennes free energy. We analyze the morphology of layers at the grain 
boundary, to better understand the mechanism of frustration between the smectic layer order and 
chirality. As the chirality increases, the layer compression energy strongly increases while the ef- 
fective layer bending rigidity is reduced due to unlocking of the layer orientation and the director. 
This results in large deviation of the layer morphology from that of Scherk's first minimal surface 
and linear stack of screw dislocations. 

PACS numbers: 61.30.Jf, 61. 72. Mm, 61.20.Ja 
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Frustration causes complex structural patterns in a va- 
riety of materials and over different lengthscales. To re- 
solve the frustration, the equilibrium patterns often con- 
tain topological defects. Liquid crystals are some of the 
richest materials providing a variety of frustrated defect 
phases. Defects are inevitably formed by frustration be- 
tween the smectic order and chirality, because the peri- 
odic layer structure and continuous helical structure are 
geometrically incompatible. As a result, the frustration 
gives rise to a set of screvif dislocations. The simplest 
stable phase containing such defects is the twist-grain- 
boundary (TGB) phase, the existence of which was pre- 
dicted using the analogy with the Abrikosov vortex lat- 
tice of type-II superconductors 0|. In this phase, smectic 
slabs of a certain length (grains) are twisted by a certain 
angle from their neighbors and separated from them by 
a narrow region (grain boundary) in which screw dislo- 
cations are aligned. 

After its theoretical prediction and experimental 
confirmation 0, study of the TGB phase was extended 
to other chiral frustrated phases and transitions between 
them. Theoretically, thermal fluctuation changes the 
transitions between smectic-A (Sm-A), TGBa, TGBc 
and cholesteric (N*) phases from second-order to first- 
order. In parallel to the vortex liquid phase of supercon- 
ductors, fluctuation induces the chiral line (A^£) phase 
with a melted defect lattice, which was theoretically pre- 
dicted and experimentally confirmed 0, I3- Many 
other chiral frustrated phases have been found and re- 
ported 

On the other hand, the spatial structure of the TGB 
phase, especially that of a grain boundary, still has much 
to be understood. Kamien and Lubensky Q showed 
that the grain boundary structure is well described by 
a linear stack of screw dislocations (LSD) if the twist an- 
gle a is close to zero. For small a, the layer structure 
is also approximated by Scherk's first minimal surface. 
However, the minimal surface defined by vanishing mean 
curvature H is achieved when the bending elasticity is 
the only contribution to the free energy. Deviation from 
the minimal surface should occur due to the other elas- 



tic effects, namely the layer compression and the twist 
Frank energies. Also, spatial variation of the density 
near the TGB defect core is not considered in the pre- 
vious analytic studies, and the interplay between these 
effects remains to be investigated 0, 0, . Study of 
the TGB structure for large twist angles should be also 
helpful for understanding the structure of more complex 
phases, such as the iV£ and smectic blue phases 0,^3- 

It might be worth mentioning that the TGB structure 
is also exhibited by other materials such as twisted lamel- 
lae of block copolymer melts and defect-containing 
Turing patterns in a reaction-diffusion system jl4j . The 
mechanism that determines the defect structure depends 
on the material, and the uniqueness of the liquid crystal 
TGB structure is not yet clarified. 

In this Rapid Communication, we study the core struc- 
ture of a grain boundary by numerical minimization of 
the Landau-de Gennes free energy. Thus, compared to 
previous molecular simulations of the TGB phase p^ll^ . 
our results allow more direct comparison with analytical 
results. Also, the previous simulations suffer from severe 
finite-size effect because the size of the simulation box 
and the grain size £b are generally incommensurate with 
each other. By focusing on a single grain boundary, we 
become free from this problem and can study the core 
structure with a higher resolution. The neglected grain 
boundary interaction would not affect the core structure 
except in the vicinity of the TGB-iV* transition. The 
difference between the simulation result and the model 
surfaces (i.e., Scherk's first surface and the LSD) is ana- 
lyzed as a function of the twist angle a. We discuss the 
reason for the deviation in terms of the coupling between 
the smectic order parameter and the director. 

The free energy we utilize is the covariant Landau-de 
Gennes model expressed as 0, IT^ 
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FIG. 1: Single twist grain boundary structure. The broken 
and gray lines indicate the parallel screw dislocations and the 
smectic layers in the grain boundary plane {z = 0). The solid 
and dotted lines are the layers a,i z —> ±00, respectively. 
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where 5* is the smectic (complex) order parameter com- 
bining the density modulation (oc Re and layer dis- 
placement, and n is the director. The double-well poten- 
tial i^D.w. with the dimensionless temperature r controls 
the order-disorder (TGB-7V* or Sm-A-A^*) transition. 
The interaction part Fint with the coupling constant B 
fixes the layer thickness d to 27r/go- The third term is 
the Frank elastic energy. Although the ratios between 
the Frank elastic constants Ki affect the dislocation core 
structure |3| , we set all the Frank constants to the same 
value K and focus on the role of chirality and coupling 
between ^I^ and n. The correlation length of the order 
parameter '5 is defined as ^ = y^J/M , the director pen- 
etration length A = Kg/ B\t\ and the Ginzburg param- 
eter K = A/^ = ^/gK /qqE. The TGB phase is stable only 
when K > l/\/2 according to the mean- field theory 1]. 

The numerical minimization is performed in a L^; x 
Ly X Lz simulation box, with the y- and z-axes iden- 
tified with the direction of dislocations and the heli- 
cal axis, respectively. To fix the transverse dimensions 
Lx and Ly, we fix the twist angle per grain bound- 
ary a as an independent parameter, instead of the chi- 
rality ko- Then, two smectic slabs sandwiching a sin- 
gle grain boundary satisfy a two-dimensional (2D) crys- 
talline symmetry with the periods £x — d/ sin(Q;/2) and 
iy = d/cos(a/2) (see Fig. 1). Thus we can assume the 
periodic boundary condition in the transverse directions. 
We set Lx = £x and Ly = £y to save computation time. 
The layer orientation changes by the angle a along the 
helical axis, which is imposed as a boundary condition 
at z = and z = as follows. The director is set to 
n — n± = (±sin(Q!/2),cos(a/2),0) sX z = L^ and z = 
respectively, while the smectic order parameter at these 
boundaries are connected by inversion with respect to 
the plane x = Lx/i: ^i{x,y,0) = '^{Lx -x,y,Lz). The 
latter is compatible with the 2D crystalline symmetry in 
the xy plane. For finite L^, these boundary conditions 
induce interaction of the grain boundary and its images. 
However, the interaction is expected to decay exponen- 
tially Jl and hence our boundary condition gives a good 
approximation if Lz/2 > ^, A. We will use Lz = 2d in 



FIG. 2: Snapshot of the grain boundary structure at r = 
-0.02,3 = 1,-B = 0.2, K = 0.02 and a = 50°. Plotted are the 
isosurface Re^' = and the director n (arrows). 



the simulation. 

The free energy is minimized by solving the time- 
dependent Ginzburg-Landau (TDGL) equations. 
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where T^J, and r„ are appropriate kinetic constants and 
the factor 1 — nn in © ensures that |np = 1. For the 
initial condition we take two flat smectic slabs with sinu- 
soidal order parameter profiles and layer normals iden- 
tical to n± for z > Lz/2 and z < Lz/2 (respectively). 
Since we fixed a instead of ko, the chirality is determined 
by minimizing Fprmik with respect to k^ as 



/co = j dr{n-y xn). 
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which is calculated from the director configuration. In 
the equilibrium, it is a function of the temperature r 
and the twist angle a: ko — f(T,a). Inverting this re- 
lation, we can find the equilibrium twist angle for given 
temperature and chirality as a = .g(r, ko). We use the 
parameter set t = -0.02, g = 1,B = 0.2, K = 0.02 and 
d = 16 with the unit mesh size Ax = 1, unless otherwise 
stated. For this choice, the correlation length and the 
penetration length are ^ = 0.44c? and A = 0.80d so that 
the conditions Lz/2 > £,,X and k > 1/V2 are satisfied. 



The kinetic coefficients are chosen as F 



0.1 



with the time step At = 1. Our criterion of equilibra- 
tion is that the free energy difference Fi{t + to) — Fi{t) 
is lower than 0.01% of the characteristic amount of Fi 
for Vi, where the index i indicates each free energy com- 
ponent and the relaxation time to is determined by the 
fitting F{t) ^ exp{—t/to) + const, at the early stage 
t < 2000. To prevent a trapping by local free energy 
minima, we also added a small noise which is gradually 
reduced to zero. The equilibrated grain boundary struc- 
ture for a — 50° is shown in Fig|21 

We first compare the obtained layer structure with 
Scherk's first surface. The deviation from the minimal 
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FIG. 3: The spatial average of the squared and dimensionless 
mean curvature (^{HdCf) versus the twist angle a at r= (x) 
-0.005, (0) -0.02 and (•) -0.05. 



surface can be estimated using the spatial average of 

1 



(8) 



This measure has more direct physical meaning than the 
previous ones 13] because (jHJ is proportional to the bend- 
ing elastic energy of layers. The mean curvature is cal- 
culated through the layer normal m as H — W ■ m, while 
m is calculated through the phase $ of the smectic order 
parameter 5* as m = V<i>/|V<i>|. In this way we can com- 
pute the mean curvature at every point, which is very 
important for averaging out the effect of mesh size. 

The twist-angle dependence of the layer curvature is 
shown in Fig|21 For small a, the TGB layer is close to 
Scherk's first surface, agreeing with the analytical calcu- 
lation However, {H'^) grows roughly linearly as a 
function of a. For large a, the mean curvature is not 
small even compared to the inverse layer thickness 1/d. 

To understand the origin of the deviation, we note that 
the layer bending elasticity has a contribution from Frank 
elasticity, through the coupling between and n. To see 
this, it would be instructive to rewrite the interaction 
term Fint into the form 



B 



dr 



|*||V$|(m - n) - I^KIV^I - qo)n 
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The first term homogenizes the order parameter ampli- 
tude, the second term locks the layer normal m and the 
director n, and the third is the layer compression term 
that adjusts the layer thickness. Thus the coupling be- 
tween and n is divided into the locking and layer com- 
pression terms. If the locking effect is dominant, m and 
n will be identical. Then the splay term in the Frank 
elastic energy is converted into the layer bending energy, 
as {K/2) j dr{V ■ nf = {K/2) J drH^. However, if the 
locking is weak compared to the twist free energy (cx kg) 
and/or the layer compression energy, an unlocking of m 
and n should occur. Then the splay term contributes 
less to the effective layer-bending energy and {H"^) can 
easily deviate from zero. 

As a increases, the factor (|V<I>| — qq)^ in the layer 
compression energy greatly increases (Fig0Jb)), as it is 
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FIG. 4: (a) Twist-angle dependence of (•) the dimen- 
sionless locking energy /; = ((;/|r|) {|*I'||V<E'|(m — n)d} 
and (0) the dimensionless layer compression energy fc = 
(ff/kl) {|*o|(|V$| -qo)d}'^. (b) The layer compression fac- 
tor fco — {cos(q/2) — 1}^ (see text). 



approximated as Qq {cos(q;/2) — 1} at the grain bound- 
ary core. Thus the frustration between layer compression 
and locking terms increases and leads to the unlocking of 
m and n. The spatially-averaged angle between them is 
over ~ 20° for the standard parameter set and a = 90°. 
Note that the amplitude | ^' | near the grain boundary de- 
creases as a increases at higher temperature, as shown in 
FigI3 However, it does not affect the ratio between the 
two free energy contributions and hence is not a major 
cause of the deviation from the minimal surface. 
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FIG. 5: The squared and dimensionless order parameter 
jv|/|^/(|r|/(;) in the cross-section y = const, at r = —0.005 and 
a— (a) 30, (b) 60 and (c) 90° . The contour lines are drawn at 
0.1, 0.2, . . ., and 1.0. For large a, a significant melting of the 
smectic order is observed even far from the dislocation cores, 
while it is not seen at lower temperature. 



We further tested the role of director in two ways: (i) 
by decoupling the director from ^ by turning off the 
Frank elasticity, and (ii) by adding an extra locking term 
D{m — n)^ to the free energy. In the unlocking limit 
(i), {H"^) increased about five times, which proves Frank 
elasticity to be the major contribution to the layer bend- 
ing rigidity. With (ii), (i?^) is reduced to one half for a 
weak extra locking (which reduces ((m — n)^) by only 
12 %). These results confirm that the director unlocking 
plays an important role in the deviation from the minimal 
surface. 

Next we look at the temperature dependence. While 
we varied the temperature r (roughly proportional to 
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FIG. 6: Deviation of the layer normal from that of (•) 
Scherk's first surface and (0) LSD. Each data set is smoothed 
by averaging over three consecutive points. 

(l^'P)) by a factor of ten, the deviation (ff^) showed 
only a small change. This also supports that the devia- 
tion is controlled mainly by the ratio between the above 
two free energy contributions. At high temperature, the 
twist free energy becomes important compared to the 
locking term, which is reduced due to the prefactor \'^\. 
The deviation for a large twist angle is larger for higher 
temperature. This can be explained by the melting of 
the smectic order in the grain boundary, not only near 
the center of the dislocation core (Fig[SJl. 

Next, we compare the obtained layer structure with 
both Scherk's first surface and the LSD. To this end, we 

use the measure (\m — mj.|^^, where rrir = V<i>r/|V$r| 

is the normal vector of the reference surface defined by 

the phase function 

-If / a\ ~| a TT 

= tan |tan yx sin —j tanh zj+y cos ^ '^i (1") 

where z — (zsinck)/2 for Scherk's first surface and 
z = 2:sin(a/2) for the LSD. The resultant deviations 



(^\m - mscherkl^^ and (\m ~ rriLSDl'^^ are plotted in 
FigEl We see that the LSD gives a better approximation 
of the TGB structure for a large twist angle, while the 
difference is negligible for small a. 

Finally let us compare our results with the TGB struc- 
ture in twisted layers of block copolymer melts 13] . In 
block copolymer melts, deviation of the intermaterial di- 
viding the surface from the minimal surface is caused 
by the packing frustration 19], which corresponds to 
the layer compression energy of liquid crystals. A self- 
consistent-field-theoretic calculation |13| shows that the 
LSD is a better model than Scherk's first surface, as in 
our case. However, the measure of deviation in [13'] has 
the dimension of the squared layer displacement, which 
causes an apparent rise of the deviation as a — > 0. To 
compare with our results, their measure must be multi- 
plied by the square of the characteristic wavenumber pro- 
portional to sin^(a/2) ^ o? . Then the deviation from the 
LSD and minimal surfaces should converge to as a — *■ 0. 

In summary, we have investigated the structure of a 
single grain boundary in the TGB phase. Sizable de- 
viations from the model surfaces are obtained and in- 
terpreted by increase of the layer compression energy 
and unlocking of the director and the layer normal. At 
higher temperature and large twist angle, the smectic or- 
der melts even far from the dislocation core. We hope 
that the core structure may be observed experimentally, 
using the same kind of technique as used in the study 
of defect cores in the smectic phase |23]- We plan to 
simulate the 7V£ phase in the future, for which the weak- 
ening of effective layer-bending elasticity may have some 
significant effect. 

We are grateful to Toshihiro Kawakatsu, Helmut 
Brand, Jun-ichi Fukuda, and Tomonari Dotera for vari- 
ous helpful comments and discussions. 
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